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ABSTRACT 

We present optical and submillimetre polarimetry data of the Bok globule CBS. We 
also present optical polarimetry and submillimetre continuum data of the Bok globule 
CB246. We use each set of polarimetry data to infer the magnetic field orientation 
in each of the clouds. The optical data can only trace the field orientation in the low 
density edge regions of clouds, because if the extinction is too high then no optical 
emission is transmitted. The submillimetre data can only trace the field orientation 
in the high density central regions of the clouds, because current submillimetre po- 
larimeters are only sensitive to high column densities. 

It has previously been found that near-infrared polarisation mapping of back- 
ground stars does not accurately trace the magnetic field in dense cloud regions, and 
hence that the grains responsible for near-infrared polarisation are under-represented 
in those regions. This may be due to a lack of aligned grains in dense regions. We test 
this by comparing the field orientations measured by our two independent methods of 
optical and submillimetre polarimetry. We find that the field orientation deduced from 
the optical data matches up well with the orientation estimated from the submillimetre 
data. We therefore claim that both methods are accurately tracing the same magnetic 
field in CB3. Hence, in this case, there must be significant numbers of aligned dust 
grains in the high density region, and they do indeed trace the magnetic field in the 
submillimetre. 

We find an offset of 40±I4 degrees between the magnetic field orientation and the 
short axis of the globule. This is consistent with the mean value of 3I±3 degrees found 
in our previous work on prestellar cores, even though CBS is a protostellar core. CB246 
is a prestellar core, and in this case the offset between the magnetic field orientation 
inferred from the optical polarisation data and the short axis of the core inferred from 
the submillimetre continuum data is 28±20 degrees. Taken together, the six prestellar 
cores that we have now studied in this way show a mean offset between magnetic field 
orientation and core short axis of ~30±3 degrees, in apparent contradiction with some 
models of magnetically dominated star formation. 

Key words: stars: formation - ISM: clouds ~ ISM: dust, extinction - polarimetry 



1 INTRODUCTION 

Many models of star formation predict that mag- 
netic fields play some role in the star formation 



process jMestel fc Spitzerl Il956l: iMouschoviasI Il976a 



Nakano fc Nakamural 1 19781: ICiolek fc Mouschoviad 1 199 



Ostriker. Gammie. fc Stonel 1 19991 : iLi fc Nakamural l2004h . 
However, measuring magnetic fields in space is notoriously 
difficult. Zeeman measurements allow one to measure the 
line-of-sight field strength using the line splitting of differ- 
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ent electronic magnetic moment states in the pres ence of 
a magnetic field ([Zeeman ,1897, : Verschuur 1969; Crutcherl 
ll999l : lTroland fc Crutcherl2008f ). Faraday rotation of a back- 
ground polarized source seen through a region containing a 
magnetic field can also be used to measure the line-of-sight 
field strength through that region (|Burnl 19661 : ISokoloff et al.l 

n"998,y 

The orientation of the magnetic field in the plane 
of the sky can be measured via the effects of aligned 
dust grains. The method of grain ali gnment has been a 



matter of debate for many years (e.g. Davis fc Greenste i nl 
Il95ll : I Jones fc Spitie3l 19671 : |Purcelllll979l : lLazarianlll995al lbl: 
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iDraine fc Weingartnej[l996l ). However, the various methods 
are generally lum ped together under the nam e of the Davis- 
Greenstein effect (jPavis fc Greensteinlll95lh . 

All methods predict that the grains align with their 
long axes perpendicular to the magnetic field. Hence in 
the optical and near-infrared, where the polarization is 
caused by preferential extinction by alig ned el ongated grains 
of a background unpolarized source iHalll fl949l : Iffiltnej 
the polarization orientation is parallel to the field. 
In the far-infrared and submillimetre, where the polar- 
ization is caused by preferential emiss ion from aligned 
elongated dust grains (.Dennison 197 7: Cudlip et al.l 1 19821 : 
iHildebrand. Dragovan. fc NovakI 19841 '). the polari zation ori- 
entat ion is perpendicular to the field orientation l|Dennisonl 
1 19771 ). 

The method of tracing the magnetic field using optical 
polarization of background starlight typically only works in 
regions of low extinction {Ay ~ 1-5). So this technique is 
mainly applicable towards the low-density edge regions of 
the clouds. This is because these data use preferential extinc- 
tion of background starlight by aligned grains to produce the 
polarisation. At higher extinctions the background starlight 
tends to be completely extingu ished and it is not possible to 
measure any polariz ation (e.g. IVrba. Strom, fc Strom|[l976l : 
IWilking et alll 19791 ) 

The method of tracing the magnetic field 
using submillimetre polarization by preferential 
emission can typically only be measured in re- 
gions of high exti n ction - Ay ~ 10-100 (e.j 
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2OOII : iKirk. Ward- Thompson, fc Crutcheij I2OO6I ). Hence 
the submillimetre data can only trace the field orientation 
towards the high density central regions of the clouds. This 
is because these data use preferential emission from aligned 
dust grains to produce the polarisation, and are only 
sensitive to very high column densities in cloud centres. 

At lower column densities, submillimetre polarimeters 
do not currently have the sensitivity to detect polarized 
emission. Furthermore, the nature of the usual data-taking 
method, using a chopping secondary mirror, ensures that 
any spatially extended components associated with the outer 
parts of the cloud are chopped away. Consequently, a com- 
bination of these two methods, of optical and submillimetre 
polarimetry, in theory provides a way of tracing the mag- 
netic field from the low extinction outer parts of molecular 
clouds, to the high extinction, compact, central, dense parts 
of clouds. 

It has previously been found d Goodman et all Il995l ) 
that near-infrared polarisation mapping of background stars 
does not accurately trace the magnetic field in dense cloud 
regions. This is based on the observation that near-infrared 
polarisation levels, caused by preferential extinction by 
aligned grains , do not appear to ris e sufflciently in denser 
cloud regions ijGoodman et al.ll 19951 ). 

This is not an effect of extinction preventing the in- 
frared emission from escaping the cloud. Instead, these au- 
thors concluded that the grains responsible for polarising 
backgro und starlight are unde r-represented in dense cloud 
regions ijGoodman et al.lll995l l. There are two possible ex- 
planations for this: either there is a lack of aligned grains in 
dense regions; or that only a small fraction of the dust grain 



size population in dense regions lies in the correct size range 
to polarise near-infrared emission from background stars. 

In this paper we attempt to test the first of these hy- 
potheses, whether or not there are significant numbers of 
aligned grains in high density regions. We do this by compar- 
ing submillimetre polarization measurements in high-density 
regions with optical polarization measurements in lower- 
density regions of the same cloud for a case in which both 
sets of measurements have been taken. 

The submillimetre data cannot be affected by lower den- 
sity foreground, or extended, field components, since they 
are not sensitive to these components. If the submillimetre 
and optical polarisations agree, then it would be difficult 
to avoid the conclusion that the magnetic field is continu- 
ous from low to high density and that the dust grains in 
the high-density regions are tracing the magnetic field in 
those regions. This would lead to the conclusion that there 
are significant numbers of aligned grains in the high-density 
regions. 



2 OBSERVATIONS 

Optical polarization measurements were made at the 1.2- 
metre telescope at Gurushikhar, Mount A bu, India, at lat i- 
tude = +24:° 36', longitude = +72° 43' (|Sen et al.l[200ol ). 
The Cassegrain focal station was used, giving an //13 
focus. The instrument used was the Imaging Polarime- 
ter, IMPOL, of the Inter University Centre for Astron- 
omy and Astrophysics (lUCAA) , Pune, India. A des crip- 
tion of the instrument is given bv lSen fc TandonI (|l994l ) and 
iRamaprakash et al.1 l|l998f ). 

The design of the IMPOL p olarimeter was based 
on the Durham Polarim eter (|Scarrott et al.l Il983l : 
IWard- Thompson et al.1 Il985l ). It uses a WoUaston prism 
and half-wave plate to observe two orthogonal polarization 
states simultaneously. This removes the effects of sky 
polarization. A CCD camera was used to image the data. A 
broad-band white light filter was used, which corresponds 
to the wavelength range 0.45-0.7 iim. The observations 
were carried out with typically a 500-sec exposure time for 
each object. The pol arimetric d ata reduction procedure is 
described in detail bvls^ et al. I lipoo). 

We here report on optical polarization data of 
Bok globules CB3 a nd C B246, from the catalogue of 
IClemens fc Barvainid ll 19881) . These data were first pub- 
lished bv ISen et al.1 (|200G| ) as part of a larger study. The 
data were taken on 1997 December 23-25. The polarization 
was measured of stars background to the clouds and seen 
through the cloud periphery. The stars are generally of vi- 
sual magnitude 12 or fainter. 

Maps of the optical polarization vectors for these glob- 
ules were mostly produced for the outer parts of the clouds, 
where the optical extinction was not so high. CB3 exhib- 
ited polarization values generally below 2% , whereas for 
CB246 it was a little higher. A detailed analysis of the na- 
turc of these polarization values was presented bv lSen et al.l 
(2005). 

The submillimetre observations were made at the James 
Clerk Maxwell Telescope (JCMT), Mauna Kea, Hawaii, 
using the Submillimetre Common User Bolometer Array 
(SCUBA). SCUBA contained a pair of hexagonal, close- 
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Right ascension 

Figure 1. Polarization map of CBS. The thin vect ors are the optical pol arisation vectors from lSen et al. I 1I2OOOI) . The thick vectors are 
the submm polarisation vectors from SCUBAPOL ijMatthews et al]|2009l '). which have been rotated through 90 degrees to illustrate the 
magnetic field orientation in the same manner as the optical polarisation vectors (see text for discussion) . The underlying greyscale image 
is the optical image of the region, taken from the digitized sky survey (Lasker 1994), obtained using the Sky View interface (McGlynn 
&; ScoUick 1994). The contours are isophotal contours of 850-/im continuum emission as mapped by SCUBA. Note the good agreement 
between the field orientations inferred independently from the optical and submm polarisation data. 



packed, bolometer arrays that could take data at 850 and 
450 microns simultan eously, using a dichroic beam-splitter 
ijHoUand et al.l Il999l ). SCUBA polarisation measurements 
were taken wi^th th e SCUBA polarimeter, SCUBAPOL 
jCreaves et al.ll2003l ). which used a rotating half-waveplate 
and fixed analyser. 

The waveplate was stepped through 16 positions (each 
offset from the last by 22.5°) and a Nyquist-sampled image 
(using a 16-point jiggle pattern) was taken at each waveplate 
position (Holland et al. 1996; Greaves et al. 2003). The ob- 
servations were carried out whilst chopping the secondary 
mirror 120 arcsec in azimuth at 7 Hz and synchronously 
detecting the signal, thus rejecting sky emission. The in- 
tegration time per point in the jiggle cycle was 1 sec, in 
each of the left and right beams of the dual-beam telescope 
chop. The total on-source integration time per complete cy- 



cle was 512 sec. Only 850-micron data have been used in 
this paper. The instrumental polarization (IP) was removed 
using the measured value of Greaves et al. (2003). This 
was 0.92±0.05% at a position angle of 163±2 degrees. This 
was measured on planets and found to be stable with time. 
The main cause of the IP was believed to be the telescope 
wind blind. SCUBA and SCUBAPOL suspended operations 
in 2006, pending the commissioning of the next generation 
SCUBA2 instrument. 

iDi Francesco et al.l (|2008l ) and [Matthews eFal] (|2009l ) 
have released near-complete legacy datasets for SCUBA and 
SCUBAPOL respectively. The SCUBA p olarisation data 
presen ted in this paper come from the [Matthews et al.l 
j2009t ) legacy data set. The data were originally taken on 
19 98 August 10. The Stok es I, Q & U maps were calculated 
bv lMatthews et al] (|2009l ) by fitting a sinusoid to the inten- 
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Table 1. Position angles of the magnetic field vectors shown in 
Figure [T] across CB3. The submm polarisation angles have been 
rotated by 90 degrees to indicate the orientation of the magnetic 
field, for direct comparison with the optical polarisation angles 
(see text for discussion). 



sity measurements from the different waveplate angles. The 
Stokes maps were then converted by us into tabulated lists 
of polarisation vectors. We binned the I, Q & U maps to a 
resolution of 20 arcsec, and then calculated the polarisation 
vector at each pixel. Finally, a signal-to-noise ratio cut was 
applied to remove spurious vectors. The remaining vectors 
are plotted in Figure 1. 

The continuum submillimetre data were taken sepa- 
rately from the polarisation data. For CB3 the data were 
downloaded from the SCUB A data archive w h ile th e data 
for CB246 were taken from iLaunhardt et al.1 l|2009t ). The 
data were reduced in the usual ma nner (see, for exam- 
ple, iNutter fc Ward-Thompsonll2007fl us ing the SURF soft- 
ware package i Jenness fc Lightfootll200ol ). The submillime- 
tre zenith opacity for atmospheric extinction removal was 
determined by compa rison with the 1.3-mm sky opacity 
l| Archibald et al.ll2002l ). 



3 RESULTS 
3.1 CB3 



CBS was catalogued by I Clemens fc BarvainisI (|l988l ). who 
noted it has an angular si ze of 6.7x5. 6 arcmin. It is 
also known as Lynds 594 (|Lvndsl 19621'). It lies at an 
appro ximate distance of 2.5kpc (|Launhardt fc Hennind 
1 19971 ). It contains at least two young stellar ob- 
j ects (YSOs), one of which is a C lass source 
l|Andre. Ward- Thompson, fc Barsonvlll993|) driving a chem - 
ically active bipolar outflow ( Codella fc Bachillerl Il999l ). 
Additionally there appear to be more than 20 very red 
K-band sources i n the field clustered around the core 
l|Launhardt et all Il998l ). This cluster was compared to 



the su bmillimetre emission bv lHuard. Weintraub. fc Sandelll 
l|2000l ). CBS is sometimes referred to as a site of intermediate 
mass star formation (|Codella fc Bachillerl [l999l :lKi ein et al.l 
2QQ0. ). It is coincident with the IRAS source 00259+5625, 
which has a luminosity of ~9S0 L© (| Codella fc Bachillerl 
Il999l) . There is also a water maser source associated with 
this globule (|Scappini. Caselli. fc Palumbd lT991). 

Figure [T] shows the optical and submillimetre polarisa- 
tion maps of CBS, superposed on a contour map of submil- 
limetre continuum emission and a grey-scale image of optical 
continuum emission. The region of brightest submillimetre 
emission coincides with a blank area in the optical image, 
as expected for a dark dust cloud. Nevertheless, background 
stars are visible through the edge of the cloud. The optical 
and submillimetre polarisation vectors mostly appear to lie 
in a similar orientation. 

We find the E(B— V) colours of the background stars to 
all lie in the range between 0.1 and 0.6 (Sen et al., in prep.). 
Taking the canonical ratio (Whittet et al., 2008, and refer- 
ences therein), we therefore estimate that these stars have 
Av O.S-1.8. In the region covered by the outermost sub- 
millimetre contour, the mean column density is ~10^^cm~^ 
(Launhardt et al., 1997). This corresponds approximately 
to Av ~ 10 (Frerking et al., 1982). This is precisely the 
level of extinction in which Goodman et al., (1995) found 
near-infrared polarisation levels lower than expected. 

Table [T] lists the positions and position angles of the 
submillimetre vectors across CBS. The angles have been ro- 
tated by 90 degrees to show the orientation of the magnetic 
field. The mean and standard deviation of the position an- 
gle of the magnetic field derived from the Q and U values at 
each of the 10 positions is 69±14 degrees. 

Table [1] also lists the six optical magnetic field vectors 
which are closest to the submillimetre vectors - those optical 
vectors within an arcmin of the centre of the submillimetre 
emission. The mean and standard deviation of the position 
angle of these six optical vectors is 77±15 degrees. Vector 
number 20 is seen to be significantly different from the other 
five vectors. This may be due to it being a foreground star. 
We note that it does have a lower than average value of E(B- 
V), consistent with this explanation (Sen et al., in prep.). 
If we exclude this vector and calculate the mean of the re- 
mainder, we obtain a mean of 72±4 degrees. 

Therefore, we see that the field orientation as estimated 
from the submillimetre polarisation is in very good agree- 
ment with that estimated from the optical polarisation. Con- 
sequently, we deduce that the magnetic field is continuous 
from the low to the high density regions, with a mean po- 
sition angle of 71±14 degrees. Furthermore, the dust grains 
in the high density region appear to be tracing the magnetic 
field well. 

The submillimetre contours on Figure [1] show that this 
Bok globule is elongated. The same elongation is seen in var- 
ious molecular line tracers of this core, such as CS (Codella 
fc Bachiller 1999). We estimate that the position angle of 
the long axis of the globule is 21±5 degrees. This is at an 
angle of 50±14 degrees to the magnetic field orientation. 
Models of magnetically dominated star formation predict 
that the magnetic field should lie along the short axis of the 
star-forming core (e.g. ICiolek fc MouschoviasI 1 1994 Il998l : 
iLi fc Shu|[l99^ . In this case we see an offset of 40±14 de- 
grees from the short axis. This is similar to results we have 
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Figure 2. Polarization map of CB246. Details as in Figure 1. No submm polarisation data were available for this region. However, note 
how the field orientation inferred from the optical polarisation data lie offset from the short axis of the submm continuum emission from 
the core. Note also how both core and field appear to twist (see text for discussion). 



seen previously in prestellar cores (|Ward-Thompson et al.l 
I2OOOI : iKirk. Ward- Thompson, fc Crutcheij|2006l ). We return 
to this in section |4] below. 



3.2 CB246 



CB246 was catalogued by I Clemens fc Barvainij ()l988h . 
who noted it has an angular si ze of 7.8x4.5 arcmin. It 
is also known as Lynds 1253 l|Lvndsl 19621). It lies at 
an a. pproximate distance of 140pc ( Launhardt fc Hennind 
It is a starless core, p ossibly prestellar in nature 
I Ward-Thompson et all |2007| ). IScappini fc Codellal (Il996l ) 
list this source as inactive, with neither a water maser 
nor molecular outflow. However, their detection of NH3 
and C2S made them postulate that the source was more 
advanced than the other sources in their survey. They 
followed up their survey with a targeted observation 



(|Codella fc Scappinil 1998|). wh i ch co nfirmed their earlier 
findings. iLee. Mvers. fc Tafallal (j200lh classified CB246 as 
having significant red excess, but no blue asymmetry in the 
CS spectra, indicating that the core has not begun to col- 
lapse. 

Figure [2] shows the optical polarisation map of CB246, 
superposed on a contour map of submillimetre continuum 
emission and a grey-scale image of optical continuum emis- 
sion. The region of brightest submillimetre emission coin- 
cides with a blank area in the optical image, as expected 
for a dark dust cloud. Nevertheless, background stars are 
visible through the edge of the cloud. No submillimetre po- 
larimetry data were available for this cloud, but we can infer 
the plane-of-sky magnetic field orientation from the optical 
polarisation vectors. 

We note that the vectors to the east of the core cen- 
tre have a systematically different mean position angle from 
those to the west. Similarly the core long axis orientation 
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Core 


Offset 


L1544 


29 ± 6 


L183 


34 ± 6 


L43 


44 ± 6 


L1498 


19 ± 12 


L1517B 


25 ± 6 


CBS 


40 ± 14 


CB246 


28 ± 20 



Table 2. Position angle offsets between the magnetic field orien- 
tation and the short axis of the core for the five prestellar cores 
previously measured, together with the two sources studied in the 
current work (see text for discussion). 

also appears to shift from east to west. Therefore, we com- 
pare the mean position angles of the magnetic field and core 
long axis separately in the eastern and western parts of the 
core. In the east, the plane-of-sky magnetic field lies at a 
position angle of 60±3 degrees. The long axis of the core 
in the east appears to lie along a position angle of roughly 
122±5 degrees. Hence the short axis is at position angle 
32±5 degrees. The approximate angular offset between the 
magnetic field and the short axis of the globule is therefore 
28±8 degrees for the eastern half of the core. 

In the west, the plane-of-sky magnetic field lies at a 
position angle of 94±5 degrees. The long axis of the core 
in the west appears to lie along a position angle of roughly 
157±5 degrees. Hence the short axis is at position angle 
67±5 degrees. The approximate angular offset between the 
magnetic field and the short axis of the globule is therefore 
27±10 degrees for the western half of the core. 

The weighted mean of these two offsets is therefore 28 
degrees. The mathematical error-bar on this value is less 
than 10 degrees. However, given that the core and field both 
curve by more than 30 degrees we do not claim an error-bar 
of less than 20 degrees for this estimate, as this is the error- 
bar obtained if one treats the whole core simultaneously. 
A detailed submillimetre polarisation map of CB246 could 
potentially confirm that the curvature of the field exactly 
mirrors the curvature of the core throughout its length. In 
that case a smaller error-bar could be ascribed to the offset 
value. Nevertheless, the value of 28 degrees is still consistent 
with prestellar cores we have observed previously (see next 
section). 



4 DISCUSSION 

The fact that we have observed a consistent magnetic field 
orientation in CBS for both the optical and submillimetre 
polarisation data shows that both are tracing the same mag- 
netic field. This implies that dust grains in the high-density 
regions of this globule are indeed tracing the magnetic field. 
Hence, there must be significant numbers of aligned grains 
in the high-density regions. Furthermore, the magnetic fields 
must be continuo us from low to high de nsity. Therefore, the 
previous finding (Goodman et al.||l995l ). that dust grains in 
high-density regions do not trace the magnetic field in near- 
infrared polarisation measurements, cannot be explained by 
a lack of aligned grains. Our findings are consistent with 
those of Whittet et al. (2008), who conclude from their data 



that there is no change in grain alignment efficiency from 
low- to high-density regions of dust clouds. 

In addition, these data allow us to comment on 
the influence of magnetic fields on the star formation 
process. Magnetically-dominated models of star forma- 
tion predict that a cloud should collapse first along 
the magnetic field lines, and then subsequently contract 
quasi-statically perpendicula r to the field orientation (e.g. 
ICiolek &: MouschoviasI Il994l ) . Hence these models predict 
that the magnetic field should lie along the short axis of 
a collapsi ng cloud core. 

War d-Thompson et al.1 (|2000l ) presented the first obser- 
vations of the magnetic field geometry in pre-stellar cores 
when they mapped L183, L1544 and L43 with the SCUBA 
polarimeter. Their maps showed relatively smooth and uni- 
form magnetic fields over the central core regions. How- 
ever, in each case they observed an offset between the 
short ajcis of the core and the magnetic field orientation. 
iKirk. Ward- Thompson, fc Crutcheil (120061 ) studied the pre- 
stellar cores L1498 and L1517B in a similar manner and also 
observed similar offsets. An attempt to explain these offsets 
in terms of tri-axial cores viewed from arbitrary angles was 
proposed by Basu (2000), but no magnetic model incorpo- 
rating this additional complication has been put forward. 

Table[2]lists the cores observed previously, together with 
those from this work, and the offset observed in each case 
between the magnetic field orientation and the short axis of 
the core. The weighted mean offset of the five previously ob- 
served cores is 31 ±3 degrees. We note that the two globules 
observed in the current study both have offsets consistent 
with this value. This is in spite of the fact that CB3 is not 
a starless core, but a protostellar core. CB246 is a starless 
core. The weighted mean of the six starless cores in Table [2] 
is ~30±3 degrees. This appears to be in contradiction with 
some models of magnetically moderated star formation. 



5 CONCLUSIONS 

We have used optical and submillimetre polarimetry of Bok 
globule CB3 to infer the magnetic field orientation. The op- 
tical data trace the field orientation in low density regions 
and the submillimetre data trace the field orientation in high 
density regions. We have tested a hypothesis that there is a 
lack of aligned dust grains in high density regions, causing 
near-infrared polarisation levels to be lower than expected, 
by comparing the field orientations measured by our two 
independent methods. 

We found that the field orientation deduced from the 
optical polarization data matches up well with the orienta- 
tion estimated from the submillimetre measurements. The 
inferred field lines join up from low to high extinction and 
produce a consistent and coherent picture of the whole mag- 
netic field in CB3. We saw that both methods accurately 
trace the same magnetic field. Hence we deduce that there 
are significant numbers of aligned grains in the high-density 
regions of CB3. Furthermore, these aligned dust grains ap- 
pear to trace the magnetic fields in those regions, as mapped 
in submillimetre polarisation. 

We have also compared polarisation data with submil- 
limetre continuum data in Bok globules CB3 and CB246. 
We have compared the field orientations with the core mor- 
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phologies and found an offset in each case between the field 
orientation and the short axis of the globule. These offsets 
are consistent with those seen in previous studies of prestel- 
lar cores, despite the fact that, of our two globules, only 
CB246 is a starless core. Taken together, the six starless 
cores that we have now studied in this way are beginning 
to show a consistent picture of an offset between magnetic 
field orientation and core short axis of ~30±3 degrees, in 
apparent contradiction with some models of magnetically 
dominated star formation. 
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